Possible origin of 60-K plateau in the YJ^CusOe+y phase diagram 
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We study a model of YBa2Cu3 0e+i/ to investigate the influence of oxygen ordering and doping 
imbalance on the critical temperature T c (y) and to elucidate a possible origin of well-known feature 
of YBCO phase diagram: the 60-K plateau. Focusing on "phase only" description of the high- 
temperature superconducting system in terms of collective variables we utilize a three-dimensional 
semi microscopic XY model with two-component vectors that involve phase variables and adjustable 
parameters representing microscopic phase stiffnesses. The model captures characteristic energy 
scales present in YBCO and allows for strong anisotropy within basal planes to simulate oxygen 
ordering. Applying spherical closure relation we have solved the phase XY model with the help of 
transfer matrix method and calculated T c for chosen system parameters. Furthermore, we investigate 
the influence of oxygen ordering and doping imbalance on the shape of YBCO phase diagram. We 
find it unlikely that oxygen ordering alone can be responsible for the existence of 60-K plateau. 
Relying on experimental data unveiling that oxygen doping of YBCO may introduce significant 
charge imbalance between Cu02 planes and other sites, we show that simultaneously the former 
are underdoped, while the latter - strongly overdoped almost in the whole region of oxygen doping 
in which YBCO is superconducting. As a result, while oxygen content is increased, this provides 
two counter acting factors, which possibly lead to rise of 60K plateau. Additionally, our result can 
provide an important contribution to understanding of experimental data supporting existence of 
multicomponent superconductivity in YBCO. 

PACS numbers: 74.20.-z, 74.72. Bk, 74.62.-c 



I. INTRODUCTION 



The YB 2 Cu306+ y compound (YBCO), as discov- 
ered in 1987 by Wu and co-workers, is the first ma- 
terial that became superconducting in boiling nitrogen 
temperature^ The material contains three copper-oxide 
layers in a unit cell: two of them are separated by Yt- 
trium atom, the third - basal plane, is surrounded by two 
Barium atoms (see, Fig. QJ. The oxygen can be intro- 
duced into the basal plane and its content can be varied 
from 6 to 7 per formula (0 < y < For y = 0, all 
Obi and Ob2 sites are empty and the system is tetrago- 
nal. As the oxygen content increases, additional atoms 
occupy Obi and Ob2 sites randomly up to a critical dop- 
ing, for which tetragonal-orthorhombic (T-O) phase tran- 
sition occurs. For higher dopings oxygen in Obi and 
Ob2 sites becomes partially ordered forming (depending 
on oxygen content) one of three different orthorhombic 
phases: ortho-I (with fragments of copper-oxide chains 
in the basal plane); ortho-II (with chains in every sec- 
ond Obi site); and ortho-Ill (with alternating one empty 
and two filled chains). Finally for y = 1 (maximum 
oxygen content), all Obi sites become occupied and all 
Ob2 - empty. The oxygen in the basal plane acts as a 
charge reservoir introducing holes into Cu02 plane cop- 
per atoms. Charge concentration can be also changed 
by substituting Yttrium atoms with Calcium, however it 



seams that holes introduced in such way tend to remain 
in Cu02 layers^ Amount of oxygen also determines elec- 
tronic state of the system. For y < 0.4 the material 
is insulating, while for y > 0.4 becomes superconduct- 
ing. The critical temperature is very sensitive to oxy- 
gen doping and the temperature-oxygen amount phase 
diagram contains two characteristic plateaus at 60K and 
90K. While the latter is now interpreted as an optimum 
doping with small overdoping region, the origin of the 
60K plateau is still not fully clear. It was argued that 
it might be explained by ordering of the oxygen atoms 
within the basal plane i 4 i 5 i 6 i 7 On the other hand, it was 
also suggested that the reason might be purely electronic: 
superconductivity is weakened at the carrier concentra- 
tion of 1/8 leading to the plateau^ It is a goal of the 
present paper to investigate the influence of oxygen or- 
dering and doping imbalance on the critical temperature 
of YBCO using a model that can accommodate strong 
anisotropy and characteristic energy scales in order to 
determine a possible origin of 60-K plateau on the phase 
diagram. 

Binding of electrons into pairs is essential in form- 
ing the superconducting state, however its remarkable 
properties-zero resistance and Meissner effect-require 
phase coherence among the pairs as well. While the 
phase order is unimportant for determining the value 
of the transition temperature T c in conventional BCS 
superconductors, in materials with low carrier density 
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Figure 1: (Color online) Crystal structure of YBa2Cu306+i/. 
Grey and open circles are copper and oxide atoms, respec- 
tively. Additionally, oxygen atoms are divided into three 
groups: from fully oxidized copper-oxide layers (O p ), apical 
(O ap ), from oxygen deficient layers (Obi, Oh). Effect on oxy- 
gen doping is also shown: a) y = — Obi and Om sites are 
empty; b) y £ (0, 1) - Oti and Om sites occupied in random 
or ordered manner, depending on actual value of y; c) y = 1 
- oxygen is ordered forming chains: all Obi sites are occupied 
and Ob2 sites are empty. 



such as high-T c oxide superconductors, phase fluctua- 
tions may have a profound influence on low temperature 
properties £ Measurements of the frequency-dependent 
conductivity, in the frequency range 100-600 GHz, show 
that phase correlations indeed persist above T c , where 
the phase dynamics is governed by the bare microscopic 
phase stiffnesses^ As a result, for underdoped cuprate 
superconductors, the conventional ordering of binding 
and phase stiffness energies appears to be reversed. Thus, 
the issue how phase correlations develop is a central prob- 
lem of high-T c superconductivity. Furthermore, it was 
directly shown that integrating the electronic degrees of 
freedom out in t-t -U-J Hubbard model, which is believed 
to correctly describe strongly interacting systems, leads 
to phase-only description of superconductivity^ 

In the present paper, we propose a semi-microscopic 
model of YBCO, which is founded on microscopic phase 
stiffnesses that set the characteristic energy scales: in- 
plane J|| and inter-plane J± couplings of Cu02 layers, 
in-plane J« coupling of basal planes and inter-plane J ± 
coupling between neighboring basal and Cu02 planes. 
Additionally, the model contains a parameter rj, which 
controls anisotropy of J,, gradually turning off basal in- 
plane coupling along b direction while r\ changes from 1 
to 0. Using our previous results^ we model values of 
in-plane phase stiffnesses as a function of oxygen amount 
to reproduce YBCO phase diagram and to elucidate the 
origin of the 60-K plateau. Our approach goes beyond 
the mean field level and is able to capture both the ef- 
fects of phase fluctuations and huge anisotropy on the 
superconducting phase transition. 

The outline of the reminder of the paper is as fol- 
lows. In Section II we construct an anisotropic three- 



dimensional XY model. Furthermore, we solve it 
in the spherical approximation with the help of the 
transfer matrix method and obtain a dependence of 
the critical temperature on model parameters: T c = 

T c fjy, Jj_, J||, J ± , Subsequently, in Section III we 

elaborate on influence of oxygen ordering on T c . We 
model values of in-plane phase stiffnesses as a function of 
oxygen amount and determine values of model param- 
eters for which the YBCO phase diagram can be re- 
produced. In Section IV, relying on experimental data 
unveiling that oxygen doping of YBCO may introduce 
significant charge imbalance between Cu02 planes and 
other oxygen sites, we show that the former are under- 
doped, while the latter - strongly overdoped almost in 
the whole region of oxygen doping in which YBCO is su- 
perconducting. Finally, in Section V we summarize the 
conclusions to be drawn from our work. 



II. MODEL 

Since, in underdoped high-temperature superconduc- 
tors, two temperature scales of short-length pairing corre- 
lations and long-range superconducting order seem to be 
well separated^ we consider the situation, in which lo- 
cal superconducting pair correlations are established and 
the relevant degrees of freedom are represented by phase 
factors < ipe (ri) < 2ir placed in a lattice with nearest 
neighbor interactions. In our notation, r, numbers lattice 
sites within £-th ab plane. The system becomes supercon- 
ducting once U(l) symmetry group governing the ipe (rj) 
factors is spontaneously broken and the non-zero value 
of (e %<pl ( Ti '\ appears signaling the long-range phase or- 
der. The Hamiltonian that we consider consists of four 
parts 

H [cp] = H\\ +H ± + H\ + H' x (1) 

containing various microscopic phase stiffnesses repre- 
senting characteristic energy scales present in YBCO 
(see, Fig. EJ : 

1. in-plane coupling Jii > within Cu02 layers: 

^ii = - J n H5Z{ cos ^ fo) ~ ^ Oj)] 

+ cos [tpu+i (rj) - V3e+i (rj)}} ; (2) 

2. inter-plane coupling J± > between neighboring 
Cu02 layer: 

H± = -J±^2^2cos[ipu{ri)-<p3t +1 (Ti)]; (3) 

t i<3 

3. in-plane coupling J,, > within basal planes, which 
can be gradually turned off along b direction by 
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Figure 2: (Color online) Structure of the YBCO supercon- 
ductor. Basal plane is in the middle of the picture with in- 
plane microscopic phase stiffness J» , which can be strongly 
anisotropic when tj parameter is small (or isotropic for r] = 1). 
Basal plane coupling with neighboring Cu02 planes is given 
by J ± . In-plane and inter-plane microscopic phase stiffnesses 
of Cu02 layers are Ju and Jj_ , respectively. 



anisotropy parameter rj G [0, 1] to simulate oxygen 
ordering: 

^i'i = ~ J 'w zZ/ZYl ■> cos ( r ») _ ^ + ^ 



t i j=-x,i 
rj cos \tpt (rj) - tp e (ri + jPj } ; 



(4) 



4. inter-plane coupling J ± > between adjacent 
Cu02 layer and basal plane: 

H' x = - J' X { C0S [VM+l fa) ~ <PU +2 ( r *)] 

I i<j 

+ cos[<y9 3f+ 2(r i )-^3£ +3 (r i )]}. (5) 

The indices i,j go from 1 to Nn being the number of 
sites in a plane, i = 1, A/j_/3, where N± denotes the 
number of layers and A" = N\\N± is the total number of 
sites. The partition function of the system reads: 



Z = 



~[dipe (rj) e 



0H[<P] 



(6) 



where (3 = l/fegT with T being the tempera- 
ture. Introducing two-dimensional vectors (r 4 ) = 
[S xt (^) , S yi (r^)] of the unit length (n) = S 2 xl (r t ) + 
S yi ( r i) = 1 defined by (rj) = [cos (r») , am <pt (rj)] , 
the Hamiltonian can be expressed in a vector form and 
the partition function written as: 



Jd 2 S e {Ti)6[Sj (rO-1] 



-PH[S] 



(7) 



where the unit length constraint (S^ (rj) = 1) is intro- 
duced by the set of Dirac-<5 functions. Unfortunately, the 
partition function in Eq. © cannot be calculated ex- 
actly. However, the model becomes solvable, while the 
rigid length constraint in Eq. Q ' 1S replaced by a weaker 
spherical closure relation^ 



6 [sj ( ri ) l] 



(8) 



Introducing different microscopic phase stiffnesses for 
Cu02 and basal planes breaks translational symmetry 
along c axis, since the inter-plane and in-plane couplings 
vary with period of 3, when moving from one plane to 
another. As a result, standard way of diagonalizing 
the Hamiltonian using three-dimensional Fourier trans- 
form of variables (S( (r^) in this case) fails, because of 
the lack of complete translational symmetry: To over- 
come this difficulty, we implement a combination of two- 
dimensional Fourier transform for in-plane vector vari- 
ables 



— ikr 



(9) 



and transfer matrix method for one-dimensional deco- 
rated structure along c-axispi^ This former operation di- 
agonalizes all terms in the Hamiltonian in Eq. (QJ with 
respect to k, leaving the dependence on layer index i un- 
changed. Therefore, the partition function can be written 
in the form: 



Z 



exp < NX + - In 

2iri I 2 



+ 00 



Jd 2 S ke 



k.C 



x exp 



(10) 



where A e £ ± (k) is an element of a square N± x N± band 
matrix, appearing as a result of non-trivial coupling 
structure along c-direction: 
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where: 



die point A = Aq of / (A): 



J|| (k) = 2J|| [cos (afcc) + cos , 
J|l (k) = 2Jii [cos (ak x ) + 77 cos (fcfc^)] 



(12) 



and Ais a Lagrange multiplier introduced by represent- 
ing the Dirac-i5 function in a spectral form S (x) = 
J-m d\/2iri exp (— Ax). The problem reduces then to 

evaluation of a determinant of the Aj§ (k) matrix (for 
technical details, we refer readers to Ref. [12l|L 

The partition function in Eq. ifTTHl can be written as: 



Z = 



+ 



dX 
2-ni 



exp[-7V/3/(A)] 



(13) 



In the thermodynamic limit N — > 00 the dominant con- 
tribution to the integral in Eq. il'M comes from the sad- 



9f (A) 



<9A 



(14) 



A=A„ 



and where /(A = Ao) becomes a free energy. In the 
spherical model, the emergence of the critical point is sig- 
naled by divergence of the order parameter susceptibility: 



G- 1 (k = 0) = 0, 



(15) 



where G _1 (k) = (Sk.^S_k.^) and (...) is the statistical 
average. Eq. (tlfill determines the value of the Lagrange 
multiplier A. Consequently, the free energy of the system 
reads: 



= A_ J_ r+ - 

; /? 3/3 



+ b dk x dk y 

- f (2ir) 2 /(ab) n 



^ ({ WW 2 [A (k)] 2 } A' (k) + 2 (0 Jl) 2 A (k) 



+y {(PJ±) 2 [A(k)] 2 } |(/?J ± ) 2 [A' (k)] 2 - [2 (/3J' ± ) 2 - A (k) A' (k)" 
where functions A (k) = /3J|| (k) — 2Aq and A (k) = /3J,, (k) — 2Ao- The Lagrange multiplier reads: 



A ° = I | 4J H + 2J l'l (! + ^) + J J- + V 8 W 
Finally, the Eq. l!14t leads to expression for the critical temperature: 

2 



J±_ + 4 J,, + 2 Jjj (1-7?) 
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(16) 



(17) 



(18) 



5 



where a = \o/fl—J\\ (£ + (), a = \ /[3-J^ (£ + 77C) and 
p (e) is a density of states of the chain (one-dimensional) 
lattice given by: 



1 



1 



? 6(1 



and 6 (x) is the unit-step function. 



(19) 



III. INFLUENCE OF OXYGEN ORDERING ON 
THE CRITICAL TEMPERATURE 

The result in Eq. l[T8|l provides us with a tool to an- 
alyze the influence of anisotropy and characteristic en- 
ergy scales present in YBCO on the critical temperature. 
However, in order to describe YBCO phase diagram, it is 
necessary to connect the model parameters with amount 
of oxygen doping. In our previous paper we have pro- 
posed a phenomenological dependence of in-plane micro- 
scopic phase stiffness of Cu02 plane on charge (hole) 
concentration, which appeared to successfully describe 
properties of superconducting homologous series 4£ 



J 



■In 



1 - ok ( s - °- i5 r 



(S) = i~W 

for 6 < 0.05 or 6 > 0.25 



for 0.05 < 5 < 0.025 



(20) 

To adapt it to the present model it is necessary to relate 
oxygen amount to charge concentration in Cu02 planes, 
since YBCO phase diagram is presented as an oxygen 
doping function of temperature. This relation was ex- 
perimentally determined by Tallon, et at. and found 
out to be roughly linear in the superconducting region 
(charge concentration changes from S = 0.05 for y = 0.4 
to S = 0.17 for y = l)pi^ Consequently, the in-plane phase 
stiffness from Eq. lf2f)j) expressed as a function of oxygen 
amount reads J11 (y) — J\\g(y), where: 



q(v) = l 1 -oh(y- ^) 2 

yyg> \ for y < 0.4. 



for 0.4 < y < 1 



(21) 



Although, it is usually stated in the literature, T-0 phase 
transition and emergence of superconductivity in YBCO 
are not simultaneous in terms of dopingili Thus, we sug- 
gest the following scenario describing the phase diagram 
of YBCO (see, Fig. 0J|: with increasing oxygen doping 
(starting from y = 0), the onset of superconductivity 
is reached (y = 0.4). The critical temperature is ris- 
ing up to a point, in which T-0 phase transition occurs 
(y = 0.5). For higher dopings, oxygen chains start to 
form. That results in interactions in basal planes be- 
coming more one-dimensional, thus phase fluctuations 
rise significantly and enough to keep the critical tem- 
perature constant. Further, while all the oxygen is or- 
dered in chains, the critical temperature roughly follows 
in-plane phase stiffness Jm dependence on doping reach- 
ing its maximum value of 93K for y = 0.95 and dropping 
slightly later. 
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Figure 3: Construction of YBCO phase diagram within pre- 
sented model: curves for T c {y, rf) for disordered (r] = 1) and 
fully ordered (77 = 0) oxygen set limits for the actual phase 
diagram. For low dopings (y < 0.5) interactions in the basal 
planes are isotropic and T c (y) dependence falls on y = 1 
curve. On the other hand, for y = 1, all added oxygen atoms 
form chains and the critical temperature reaches the other 
limit of r) = 0. Plateau region is a crossover between those 
two regimes in which variable value of y keeps T c conA. J. 
Leggett 

Figure 4: Affiliation: Department of Physics, University of 
Illinois at Urbana-Champaign, Urbana, Illinois 61801, USAs- 
tant. 



In order to realize this scenario within the presented 
model it is necessary to fix values of model parameters 
J||, Jm, Jj_, Jj_ and find a dependence of ?/ on doping 
y, which would results in constant critical temperature 
in the plateau region of the phase diagram. We assume 
that doping dependence of J|| and J,, are given by Eq. 

(42J, and, for simplicity, J± = J ± . Experimental data 
for anisotropy of penetration depth in YBCO provide us 



with ratio of Jj_/J\ 
values of J11 and J, 



Kb/K 



lOOpia Furthermore, 



can be deduced from the ratio of 



expressions for the critical temperatures (see, Fig. for 
plateau (60K) and optimum doping (93 A'): 



J\\, y = 0.5, 77 = 1 



Tjjn J' y = 0.95, 77 = 



60 
93' 



(22) 



Finally, dependence of oxygen ordering on doping r\ (6) 
can be found numerically for calculated values of model 
parameters (J|| = 19.32meV, J,, = 33.81mcV, J± = J ± = 
0.58meV) and it can by approximated by the expression: 



t](S) 



(2 - 25) 15 - 2 + (2 - 28) + (2 - 25) 



7.9 



n3.6 



(23) 
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Figure 5: (Color online) Comparison of experimental phase 
diagrams of YBCO with results of the presented model for a) 
linear and b) non-linear (mo re exact) approximation of data 
of Tallon, et al. (see, Ref lT^I . 



Figure 6: Comparison of phase diagrams for various ratios of 
J|l / J|| . Plots are normalized to have the critical temperature 
of the plateau region equal to 60K. 



The resulting phase diagram is presented in Fig. (thick 
solid line) along with experimental results, for compar- 
ison. It is clearly seen that the procedure reproduces 
characteristic features of YBCO phase diagram: 60-K 
plateau, 93K maximum critical temperature with small 
overdoped region. However, trying to fit experimental 
data better, we have performed more exact (non-linear) 
approximation of data from Ref. 0] that lead to dop- 
ing dependence of in-plane phase stiffness J» (y). Cor- 
responding phase diagram is also presented in Fig. 0as 
thick dashed line. It does not provide any new quality 
increasing the width of 60-K plateau only very slightly, 
thus it is reasonable to use simpler form of Jy (y) depen- 
dence, as in Eq. (1211 . It can be also noticed that the 
model predicts 60-K plateau to be a little bit narrower 
and moved toward lower doping region than experimen- 
tal results show. It can be argued that since it is very 
hard to control the amount of oxygen in YBCO precisely, 
some measurements, especially based on multigrain pow- 
der samples, can be inaccurate. However newer results 
are based on un-twinned single crystals and seem to be 
reliable, also because of good consistence among different 
research groups. 

Unfortunately, the most serious flaw in the presented 
scenario are the actual values of model parameters re- 
quired to reproduce experimental phase diagram. It 
is necessary for basal plane in-plane microscopic phase 
stiffness J,, to be almost two times bigger than Cu02 
planes coupling Jii . For more reasonable values Jii > 
Jii, Jj_, J ±) the influence of anisotropy in basal planes 



is almost negligible. For example for J| = J« = J, 
j' ± = 0.1J and J± = 0.01J with J = 19.3mcV, the ra- 



tio of critical temperatures for systems with anisotropic 
(rj = 0) and isotropic (rj = 1) interactions in basal planes 
reads: 



T c (y = 0) 



0.949. 



(24) 



= J, 



Similarly, for J\\ = = J, J ± = 0.1J and J| 

J,', = J' ± = 0.5J and J±_ = 0.1J the ratio is 0.983. The 
influence of basal plane anisotropy is too small to notice- 
ably change the critical temperature until in-plane phase 
stiffness J,, is increased to be higher than J||. This, how- 
ever, does not seem to be reasonable from physical point 
of view. Evolution of the phase diagram with changing 
ratio of J\\/J\\ is presented in Fig. E3 When the ratio is 
being decreased, the plateau region is disappearing. To 
summarize, we find it rather unlikely that oxygen order- 
ing into chains alone can explain the existence of 60-K 
plateau. 



IV. INFLUENCE OF DOPING IMBALANCE ON 
THE CRITICAL TEMPERATURE 

YBCO can be doped not only by adding the oxy- 
gen atoms, but also by substituting three- valent Yttrium 
atoms with two- valent Calciums. This introduces holes 
directly to Cu02 planes leaving apical and chain sites 
untouched^ even though the charge concentration within 
copper-oxide planes is high enough, lack of doping of api- 
cal sites results in interplane coupling being small and the 
onset of superconductivity is not reached. This suggests 
that charge concentration within apical sites along with 
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Figure 7: (Color online) YBa2Cu307 structure with copper- 
oxide chains fully formed. 



basal plane sites and CuC>2 layers can be significantly 
different. Similarly, as Yttrium substitution can lead to 
introduction of holes exclusively into CuC>2 planes, one 
can imagine that increasing of oxygen amounts changes 
charge concentration primarily in the chain regions and 
only some of the charges are transported to the copper- 
oxide layers. In fact, this scenario has been confirmed 
by means of site-specific X-ray absorption spectroscopy? 
Because overdoping of some regions of YBCO could pro- 
vide a factor decreasing the critical temperature with 
increasing oxygen doping and thus leading to the 60-K 
plateau, we want to investigate such a possibility within 
presented model. 

First, we want to emphasize the high significance of 
apical site doping. Analyzing a projection of the YBCO 
structure on a — c plane (presented in Fig. it can be 
noticed that a set of chains differs from full copper-oxide 
layer only by a lack of oxygens between Cu02 planes. 
Thus one can expect that the microscopic phase stiff- 
nesses J ± related to apical sites and the one within basal 
planes Ju are of the order of CuCb in-plane phase stiffness 
J|| . This suggests it is reasonable to investigate the role 
of doping of apical and basal sites on the phase diagram 
of YBCO. To investigate this , we want to determine 
values of J,, and J ± that would result in T c (y) depen- 
dence observed experimentally. For clarity, we assume 
that oxygen is fully ordered in chains to study the effect 
of charge imbalance alone. It is also necessary to model 
an influence of the number of vacancies in basal planes for 
y < 1 on the in-plane phase stiffness J,, . This relation is 
unfortunately not obvious and it is hard to find any hint 
about it specific form. Thus, we start with a simple linear 
dependence h (y) =y,so: j[ effective (y) = y 1 J\\ (y) with 
7 = 1.— Later, we choose different value of 7 = 3/2 and 
notice that although specific values of J« and J ± change, 
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Figure 8: (Color online) Calculated dependence of J« and 

J x on oxygen amount for various factors due to number of 
oxygen vacancies in copper-oxide chains. Dependences are 
normalized to have maximum value equal to 1. Dashed lines 
are guides for eye extending dependences beyond region where 
calculation was possible. 



the qualitative results are the same and in reasonable 
agreement with experimental data. 

Our procedure is as follows: we assume that anisotropy 
ratio between in-plane and inter-plane coupling among 
Cu02 planes is equal to J±/J\\ = 0.0136^ r\ parame- 
ter is equal to (oxygen in basal planes is ordered along 
a direction) and values of J« and J ± are equal and of 

the order of J|| (however J,, is modified by chosen factor 
h(y) = y 7 ). Providing values of Jm, J± and 77 into Eq. 
II18J1 along with experimentally obtained critical temper- 
atures based on resistivity measurements from Ref. 0, 
we determine doping dependence of basal and apical mi- 
croscopic phase stiffnesses. The results are presented in 
Fig. 03 Using expression in Eq. (j2f)l) it is also possible 
to calculate doping dependence of corresponding charge 
concentration (see, Fig. EJ. As it is apparent, almost 
in the whole region of oxygen doping in which YBCO is 
superconducting, basal planes along with apical sites are 
overdoped while simultaneously Cu02 planes are under- 
doped. While oxygen content is increased, this provides 
two counteracting factors, which may lead to rise of 60K 
plateau. Although specific relation of J,i and J x on y is 
dependent on assumed factor h(y) = y 7 , the results are 
qualitatively similar and in reasonable agreement with 
experimental data (see stars, in Fig. EJ. 
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Figure 9: (Color online) Charge concentration within basal 
planes, apical sites and in Cu02 planes. Stars denote experi- 
mental values of hole doping in basal planes and apical sites 
(from Ref. Q). Dashed lines are guides for eye extending 
dependences beyond region where calculation was possible. 



V. SUMMARY AND CONCLUSIONS 

We have considered a model of YBa2Cu30g+ y to study 
the influence of oxygen ordering and doping imbalance on 
the critical temperature T c (y) and to elucidate a possi- 
ble origin of well-known feature of YBCO phase diagram: 
the 60-K plateau. Motivated by the experimental evi- 
dence that the ordering of the phase degrees of freedom 
is responsible for the emergence of the superconducting 
state with long-range order, we focus on the "phase only" 
description of the high-temperature superconducting sys- 
tem. In our approach, the vanishing of the superconduc- 
tivity with underdoping can be understood by the reduc- 
tion of the in-plane phase stiffnesses and can be linked to 
the manifestation of Mott physics (no double occupancy 
due to the large Coulomb on-site repulsion) that leads 
to the loss of long-range phase coherence while moving 
toward half-filled limit (y = 0). For doping y = 0, the 
fixed electron number implies large fluctuations in the 
conjugate phase variable, which naturally translates into 
the reduction of the microscopic in-plane phase stiffnesses 
and destruction of the superconducting long-range order 
in this limit. In the opposite region of large y, the onset 
of a pair-breaking effect (at the pseudogap temperature 
T*) can deplete the microscopic phase stiffnesses, thus 
reducing the critical temperature. 

To this end, we have utilized a three-dimensional semi 
microscopic XY model with two-component vectors that 
involve phase variables and adjustable parameters repre- 
senting microscopic phase stiffnesses and an anisotropy 



parameter. The model fully implements complicated en- 
ergy scales present in YBCO also allowing for strong 
anisotropy within basal planes in order to simulate oxy- 
gen ordering. Applying spherical closure relation we have 
solved the phase XY model with the help of transfer 
matrix method and calculated T c for chosen system pa- 
rameters. Furthermore, by making a physically justified 
assumption regarding the doping dependence of the mi- 
croscopic phase stiffnesses we are able to recreate the 
phase diagram of YBCO and investigate the influence 
of oxygen ordering and doping imbalance on its shape. 
We determine that characteristic features of the phase 
diagram, i.e. 60-K plateau, can be recreated by effects 
of oxygen ordering. However, the specific values of the 
model parameters for which this result is obtained seem 
to be a little bit hard to justify. Furthermore, relying 
on experimental data unveiling that oxygen doping of 
YBCO may introduce significant charge imbalance be- 
tween Cu02 planes and other oxygen sites, we show that 
the former are underdoped, while the latter - strongly 
overdoped almost in the whole region of oxygen doping in 
which YBCO is superconducting. Increasing of the oxy- 
gen content provides then a natural mechanism of two 
counter acting factors that leads to emergence of 60-K 
plateau. Additionally, our result can provide an impor- 
tant contribution to solve a controversy of the symme- 
try of YBCO order parameter, for which various experi- 
ments give contradictory answers suggesting d or s-wave 
(although the fact that YBCO is orthorhombic should 
also lead to order parameter being a mixture of s and 
d-wave)^ Furthermore, measurements of the complex 
conductivity of high quality YBCO crystals show a third 
peak in the normal conductivity at 80K along with en- 
hanced pair conductivity below ~ 60Ki22 Authors show 
that a single d-wave order parameter is insufficient to de- 
scribe the data and successfully consider two-component 
model of superconductivity in YBCO. They claim that 
it would be tempting to assign the two superconducting 
components with the associated condensates residing on 
Cu02 planes and chains, respectively, however they do 
not see any justification of such situation. Thus, our re- 
sult can provide a natural answer for plausibility of such 
a scenario, in which s-wave component comes from over- 
doped chains and d-wave one - from underdoped Cu02 
planes. 
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